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Abstract. The early Dirac proposal to model the electron as a charged membrane is 
reviewed. A rigidity term, instead of the natural membrane tension, involving linearly 
the extrinsic curvature of the worldvolume swept out by the membrane is considered 
in the action modelling the bubble in the presence of an electromagnetic field. We 
set up this model as a genuine second-order derivative theory by considering a non- 
trivial boundary term which plays a relevant part in our formulation. The Lagrangian 
in question is linear in the bubble acceleration and by means of the Ostrogradski- 
Hamiltonian approach we observed that the theory comprises the management of both 
flrst- and second-class constraints. We show thus that our second-order approach is 
robust allowing for a proper quantization. We found an effective quantum potential 
which permits to compute bounded states for the system. We comment on the 
possibility of describing brane world universes by invoking this kind of second-order 
correction terms. 



PACS numbers: 04.20.Cv. ll.10.Ef, 31.30.jy, 46.70.Hg 
1. Introduction 

The pioneering Dirac proposal [1] (and even before by Less j2]) to model the electron as 
a charged bubble, started an exhaustive study of geometrical theories of surfaces moving 
in a spacetime modelling some relevant physical systems. This spinless Dirac geometrical 
theory describes a dynamic spherical membrane in the presence of an electromagnetic 
field where the non-electromagnetic forces are described by a constant surface tension. 
Unfortunately, as it was originally formulated, the model results limited in its predictive 
power. Over the years the model has been improved by taking into account the inclusion 
of second-order correction terms built from the extrinsic curvature of the worldvolume 
swept out by the bubble [3j Hj. Extrinsic curvature terms appear in several effective 
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actions aimed to describe surfaces in diverse contexts which accommodate relativistic 
extended objects as notable realizations of interesting physical system^f] [5j [6]. In 
addition, it is well known that the inclusion of extrinsic curvature terms into geometrical 
models is reflected on the spin content of a given theory. In that sense, an effective model 
for the electron viewed as a charged spherically membrane and based on the inclusion 
of the extrinsic curvature of the worldvolume swept out by the bubble was developed by 
Onder and Tucker (OT), in contrast to the original Dirac's model [31 Hj. Regrettably, 
the quantum approach performed for this model has the drawback of an unavoidable 
transcendental system which conceals the true dynamics of the system, and complicates 
its quantization. This model was also addressed in an harmonic estimate by studying 
radial modes [I]. The present status of this approach undoubtedly needs refmements 
since it only comprises a semi-classical approximation inherited from the transcendental 
equations arising for the aforementioned cumbersome phase space treatment. We are 
aware that this model does not pretend to be a realistic model for the electron, but, 
despite its plainness, it provides sufficient complexity in certain particle physics models, 
and also it introduces interesting resemblances to cosmological brane models which 
deserves a careful analysis |S] . 

In this paper we bring back to life the model proposed in [31 Hj, retaining the 
basic physical idea of studying an elementary particle viewed as a (2 + 1) -dimensional 
gravitating electrically charged membrane in the presence of an electromagnetic field. 
Instead of the natural tension of the membrane, to balance the Coulomb repulsion we 
consider a quite subtle surface stress derived from a correction term involving a linear 
extrinsic curvature term of the membrane. This type of terms has drawn attention 
for a long time, specially in reference to the approach in the hypersurfaces context 
by Chen [9], and by Barrabes et al as a correction term to study the dynamics of 
domain walls [TU]. Also, in order to achieve our goal, we consider an inspired alternative 
Hamiltonian approach by considering non-dynamic boundary terms. We overcome the 
usual quantization shortcomings for second-order derivative systems with the help of a 
canonical transformation and then we proceed to quantize this rigid membrane model. 
It is becoming more and more important to extended the powerful techniques of the 
Dirac treatment for constrained systems in describing physics with correctness, so that 
we would like to improve the former quantum OT approach by means of a non-standard 
strategy for the classical canonical scheme relying on the Ostrogradski-Hamiltonian 
formalism. In this way, the appropriate classical phase space of the system is identified 
and worked out by the above mentioned canonical transformation which reveals the true 
nature of the system: we are left with two first class and two second-class constraints, 
the former being of physical significance while the latter are treated as identities by 
imposing the Dirac bracket. At the quantum level, a Schrodinger-like equation is found 
and the effective potential is explicitly identified. Also, we have gained enough control 

X As a matter of fact, this is the criterion for rigidity that we follow in this work: It is a correction to 
Dirac-Nambu-Goto objects in order to penalize singularities in their evolution thus favoring a variety 
of configurations much rich in geometrical and physical structure. 
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of the global properties of tiris potential. Remarkably, we are able to show that the 
model under consideration leads to bounded states. 

The paper is organized as follows. Section [2] introduces the geometrical model by 
considering an action governing the dynamics of a membrane including a second-order 
correction term via the extrinsic curvature. In Section [3j we specialize the model to an 
spherical shell, and discuss the nature of the Lagrangian components. Section [4] and 
Section [5] are devoted to the analysis of the Ostrogradski-Hamiltonian formalism for 
our system, and to the classification of the emerging constraints. In Section [6] we go to 
the reduced phase space and perform a canonical transformation which reveals the true 
nature of the classical constraints. We develop the quantum counterpart for our model 
in Section [7j Finally, we present a summary and some concluding remarks in Section [8j 



2. The geometrical model 

Consider a 2-dimensional surface, S, evolving in a Minkowski 4-dimensional background 
spacetime with metric rj^, described by the embedding x^ = X^(£ a ) where x ii are local 
coordinates for the background spacetime (//, v — 0, 1, 2, 3), £ a are local coordinates for 
the worldvolume, m, swept out by the surface (a, b = 0, 1, 2) and X^ are the embedding 
functions for m. We consider the following effective action underlying the dynamics of 
the surface S 

S[X»} = I (-a^K + Pfe^A,) , (1) 

J m 

where K is the mean extrinsic curvature of the worldvolume constructed with the 
extrinsic curvature tensor K ab = —r\ lLV n il 'D a e v h and g denotes the determinant of the 
induced worldvolume metric g a b = V^^a^b, where e M a = X M ia are the tangent 
vectors to the worldvolume; n M is the spacelike unit normal vector to the worldvolume. 
Furthermore, D a = e^aD^, where _D M is the background covariant derivative. The 
factor a is a constant related to the rigidity parameter of the surface S, and /3 is 
the form factor of the model. Further, A^(x) is the gauge field living in the ambient 
spacetime, and j a is a fixed electric charge current density continuously distributed over 
the worldvolume, responsible for the minimal coupling between the charged surface 
and the electromagnetic field [TT]. The action functional ([TJ is invariant under 
reparametrizations of the worldvolume m. Here, we regard j a as a function depending 
only on the worldvolume coordinates and, not being a true vector on the worldvolume 
it is locally conserved on m, d a j a = 0. From a more generic perspective, the charge 
current density may be thought of a new dynamical variable, derivable from a suitable 
internal potential living on m, as developed by Davidson and Guendelman in [12]. We 
will not pursue this case in the rest of the paper. 

The variation of the action functional with respect to the embedding functions [13J 
leads us to the equation of motion 

an=^Lfn»e\F, u , (2) 
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which is a sort of Lorentz-force law. Here 1Z is the worldvolume Gaussian curvature, and 
F^ v = 2du l A u ] is the strength tensor of the electromagnetic field. It is worthy notice that 
j a must remain unchanged under arbitrary deformations of the surface, X — > X + 5X. 
Despite the effective Lagrangian density in equation ([T]) is of second-order in derivatives, 
the resulting equations of motion are of second-order. This is a remarkable situation 
since then we will not have to struggle with ghost-like issues in the quantum formulation 
of the model. 



3. A moving electrically charged bubble 

We turn now to specialize the defmitions of the previous section to the description of a 
spherical membrane S. From now on, we consider a background Minkowski spacetime 
described by ds 2 = —dt 2 + dr 2 + r 2 d6 2 + r 2 sin 2 9d(p 2 . We consider the membrane S to 
be orientable and topologically identical to S = S' 2 x K, that is, a bubble described by 
the following parametrization 

x» = X' J >(T,9,<f)) = (t(r),r(T),9,<f>) , (3) 
so that the induced metric on the worldvolume is ds 2 = g a bdC, a d^ b — —N 2 dr 2 + r 2 d9 2 + 
r 2 sin 2 6d(j) 2 where N = \/i 2 — r 2 , for simplicity, and the dot stands for derivative with 
respect to the parameter r. It is worth mentioning that N corresponds to the lapse 
function in the ADM Hamiltonian approach for branes [Hj. The normal vector to the 
worldvolume is implicitly defined by g^ u n^e u a = 0, and g^ u n tJ, n 1 ' = 1, and explicitly 
reads 

n^ = l (f, i, 0,0) . (4) 

Furthermore, for this parametrization we also have the unit timelike normal vector 
i] a = -^(1,0,0) such that g a b^ a A^t = and gabV°">l b = ~ 1 where {e%,?? a } is the S basis 
viewed from S into m and A, B = 1,2, (see [H] for more details). The non-vanishing 
components of the extrinsic curvature tensor take the following form 

K\=^(ri-ri), K\ = K% = -^. (5) 

These components are useful to compute the associated mean extrinsic curvature and 
the Gaussian curvature given b}||] 

K = — - [r(rt - f i) + 2iN 2 } , (6) 

2/ 

n = tfv? - rt) + iN 2 } . (7) 

By invoking the electromagnetic potential on the spherical shell to take the specific 
form A^ = (— -, 0, 0, 0) where q is the total electric charge on the shell, and fbdng the 
electric current as j T = q sm9, the equation of motion ^ reduces to 
d fr\ N 2 ( p N 2 f3q 2 \ 

dr\ij 2ri 3 \ 2ar 2 J ' { ' 

§ The intrinsic curvature can be computed from the contracted Gauss-Codazzi integrability condition, 
K = K 2 - K nh K ab . 
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From the action ([T| we see that the effective Lagrangian density specialized to the 
bubble reads 

c = sine {~w ^ ^ - ^ + ^ } + sL jTiAo - (9) 

Thus, our effective model for the electron, in terms of an arbitrary parameter r, is 
reduced to 

S = 4vr J drL{r,f,f,i,i), (fO) 
where the Lagrangian is given by 

L = -ajp - rt) ~ 2ar i ~ l 3 ^- ( n ) 

In addition to the velocities t and r, this Lagrangian depends also on their corresponding 
accelerations t and f. So we are dealing with a genuine second order derivative theory. 



Note that (II ) can be rewritten as L = L& + L^, where 

d 



dr 



a r arctanh 



(12) 



j j —2ari — (3 — , (13) 

are respectively a boundary term and a true dynamic term. As customary, the boundary 
term can be neglected without affecting the bubble evolution in time. This issue was 
addressed in jl] by taking into account a semi-classical harmonic oscillator approach 
due to the high degree of difficulty in the bubble evolution encountered. However, our 
treatment will rely on considering explicitly both terms, the boundary and the dynamic, 
confronting us with a Lagrangian depending up to the accelerations, hence evoking an 
Ostrogradski-Hamiltonian formalism for the treatment of the model. 

4. Ostrogradski-Hamiltonian approach 

The more general aspects for a Hamiltonian analysis of an extended object have been 
presented in detail in [H]. Here we will apply those results. 

The highest conjugate momenta to the velocities {i,r\ are given by 

P ' = 1 X=*W < 14) 
* = £ = < 15) 

such that the highest momentum spacetime vector can be rewritten as 

Note that the momentum P M is directed normal to the worldvolume. This is a general 
issue for this type of brane models 
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The conjugate momenta to the position variables {t, r} are 
8L d fdL\ 2ari 2 a 2 

dL d f 9L\ 2arir 
Pr= & : ~^\d?) = ^ 2 ~ 1 [ ' 

respectively. Important to note is the fact that both momenta, p t and J) r , £1X6 from a 
totally different nature. Indeed, while the momentum p t is not influenced at all by the 
surface terms the momentum p r is obtained by two contributions: p r coming from the 
ordinary dynamical theory (L^) and p r coming from the boundary term (Lf,). In this 
way, we can denote the momentum p r as [E] 

p r :=p r + p r , (19) 



where 



and 



lartr , ( f \ 

N2 + 2ar arctanh i j ) , (20) 



— 2ararctanh ( - ) . (21) 



It is crucial to recognize (20) as the canonical momentum worked out in [3], while (21) 
stands for the momentum conjugated to the r(r)-variable when considering as the 
Lagrangian only the surface term. 

In order to see the geometrical nature of the physical momentum it will be 
convenient to write the kinetic momentum, ir^ = p^ — (3 q A^, as follows 

2art ■ , , 

t^-Jpr**- (22) 



The constant quantity fi, previously introduced in relation (17), is nothing but the 
conserved energy as we will see below. We can rearrange the energy expression (17) in 
order to get a master evolution equation 

N 2 + r 2 = ^N 2 r 2 H 2 , (23) 

where 7 = 7(r) = (Qr — (3q 2 ) /r 4 , and we have introduced H 2 = l/2a. This relation 
represents for the bubble model an analogous master equation emerging in the context 
of quantum geodetic brane gravity [El [El [TTJ, H8] . 

Hitherto, the appropriate phase space of the system, T = {t,r,i,f;p t ,p r ,P t ,P r }, 
has been explicitly identified. In order to complete the Ostrogradski-Hamiltonian 
programme in phase space T, we will consider the canonical Hamiltonian 

H = p r r + p t i + P r r + P t i - L 

= n r r + n t i + J K , (24) 

where we have defined 

J K := V2^N(Qr-(3q 2 ) 1/2 . (25) 
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As expected, the canonical Hamiltonian results a function only of the physical 
momentum 7r M . It may look like an unnecessary complication to write both the physical 
momentum and Hq in terms of Q, but this quantity results important for the reason 
that it is merely the conserved internal energy. In order to demonstrate this, we only 



have to consider the Poisson bracket (defined below in relation (33)) of the momentum 
Pt with the canonical Hamiltonian Hq. 

4-1- Classical equilibrium configuration 

We now briefly describe the equilibrium conditions of the system. The radius of the 



static solution, r , is obtained if we put f = f = in Eqs. fl8ft and (17). Thus, we have 



/V / s 

= 2a - (26 a) 



0Q 2 

Q = 2ar + — , (266) 
r 



where (266) is the static energy in equilibrium. It is worth mentioning that these 



conditions are gauge independent. Further, note that dVL/dr Q = implies Eq. (26o), 
which is a consequence of the minimum for the classical potential function. The 
condition (266) specialized to the electron properties and Vt = m e , lead to the resuhjjj] 



r = — and a = — t, (27) 



— and a = — r 
m e Ar^ 

where A := e 2 = 1/137 is the structure constant, e 2 := q 2 and (3 = 1/2. 
5. First- and second-class constraints 

As already mentioned, we are dealing with a second-order derivative theory because, in 



addition to the velocities, the Lagrangian (11) depends also on the accelerations. From 



the defmition of the momenta (16) we can get the primary constraints 

d = P ■ X ps , (28) 
C 2 = N P ■ n + ar 2 w , (29) 

where the dot symbol denotes contraction with the Minkowski metric|5] Also, the ~ 
symbol stands for weak equality in the Dirac approach for constrained systems [191 EOj . 
The previous constraints are supported by the completeness relation 

= itn v - rfrf + h AB e» A e v B , (30) 

satisfied by the E basis, {e' 1 a, Y', n^}, where Hab = V^^a^b is the spatial metric on 
E and e' 1 a the corresponding tangent vectors to E. The vector r]^ denotes a timelike 
unit normal vector to E [1] 



| We will work in natural units, h = c = 1. 

% Note that, in general, n' 1 is a function of the derivatives with respect to the parameter r of the 
embedding functions X. 
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The total Hamiltonian that generates time evolution is given by 

H T = H + Xid + A 2 C 2 , (31) 

where Ai and A 2 are Lagrange multipliers enforcing the primary constraints C\ and C 2 . 
Time evolution for any canonical variable z G T is thus dictated by means of 

z = {z,H T }, (32) 

where the corresponding Poisson brackets (PB) for any two functions F (z), G(z) G T is 

appropriately defined as 

rr ,_, OFdG OFdG OFdG OFdG /rn 
{F ' G}:= ^ + ^ + ^ + ^- (F ^ G) - (33) 

Note that {Ci,C 2 } = and in consequence the primary constraints ( |28| and (29) are 

in involution. As primary constraints must be preserved in time, that is, C\ ~ and 

C 2 ps 0, we are lead to the secondary constraints 

C 3 = H w , (34) 
C 4 = JVtt • n w . (35) 

The vanishing of the canonical Hamiltonian is expected as for any reparametrization 
invariant theory. Geometrically, the canonical Hamiltonian Hq generates diffeomor- 
phisms normal to the worldvolume. The secondary constraint (35) is characteristic of 
every brane model linear in accelerations [15] . The process of generation of further con- 
straints is stopped at this stage since C3 is automatically preserved under time evolution 
and the requirement of being stationary for C4 only determines a restriction on one of 
the Lagrange multipliers, namely, A 2 = — (i 2 — j /Iri. Thus, we are dealing with 
an entirely constrained theory with two primary and two secondary constraints. 

Following Dirac's programme, the set of constraints should be separated into subsets 
of first- and second-class constraints [TU [20] (see also [UI ELT] in the context of geodetic 
brane gravity). For our system we have two first-class constraints and two second-class 
constraints. We judiciously choose them as 

T\ := Ci , (36) 



i 1 



/3q 2 N< 



T2 := - V 2 "f 2 7 C 2 + C 3 , (37) 

51 := C 2 , (38) 

5 2 := C 4 , (39) 

where the J-'s and the 5's stand respectively as the first- and the second-class 
constraints. In order to eliminate the extra degrees of freedom in our canonical approach 
we must replace the PB with the Dirac brackets (DB) defined by 

{F, G}* := {F, G} - {F, S,} {S v G} , (40) 

ftf 

by Sij := {Si, S j} , (i, j = 1, 2.). We find straightforwardly in our case the matrbc 

) 1 
1 



where <S,„- stands for the inverse elements of the second-class constraints matrbc defined 
{Si,Sj} , (i,j = l 

(S tj )^4art\ ° \ ). (Ml) 
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We may therefore construct the DB, and we find 

{F, G}* = {F, G} + -r— j ({F, S\} { 
Aart 

Hence, we consider the second-class constraints to vanish strongly which helps to 



{F, G}* = {F, G} + —'-r ({F, S 1 } {<S 2 , G} - {F, S 2 } {S^ G}) (42) 
Aart 



eliminate the part proportional to C 2 in (37) leading thus to a more suitable expression 



form for T 2 . The final first-class Hamiltonian for our model is H = a A pA, (A, B = 
1,2). We turn now to the counting of physical degrees of freedom (dof). According to 
the recipe developed in [2U], the model has (8 — 2x2 — 2) /2 = 1 dof. Note that as 
we have two linear independent first-class constraints, we will have the presence of two 
gauge transformations for this brane model. We will analyze the gauge-fbdng in the 
next section. 

6. Gauge— fixing 

According to the conventional Dirac scheme, in order to extract the physical meaningful 
phase space for a constrained system we need a gauge-fbdng prescription which entails 
the introduction of extra constraints, avoiding in this way the gauge freedom generated 



by first-class constraints (36) and (37). To achieve this we will consider the following 
gauge condition 



N — 1 — Vi 2 -r 2 -1«0, (43) 



and the generalized evolution equation ( 23 ) 

(p 2 = N 2 + f 2 - 7 N 2 H 2 r 2 w . (44) 

From the geometric point of view, this set of gauge conditions is good enough since 
the matrbc ({J 7 , ^1,2}) is non-degenerate in the constraint surface. Indeed, under the 



Poisson bracket structure (33), it is straightforward to show that gauges <p\ and (p% form 



a second-class algebra with the constraints T\ and T 2 

{tpuFi} = (pi + 1, 

=0, 
{^2,^1} =2v?2, 

{^2,^2} = " A L f. 

r 

Consequently, velocities t and r must be discard as dynamical degrees of freedom. 

Next, in order to get control over the model, we implement the following canonical 
transformation to a new set of phase space variables 

N := Vi 2 - r 2 , (45) 

n^v := i(P ■ X) , (46) 

v '.— — [N (P -n) + ar 2 ] , (47) 

U u := arctanh ( - J , (48) 
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together with the transformation 



P M := P M + K, *4 n, . (49) 

Important to note is that in this canonical transformation the coordinates X fl remain 
unaltered, while the dynamical momentum p r is distinguished as the relevant momentum 
of the model. Such transformation can be physically interpreted as a Lorentz rotation 
in phase space which, straightforwardly, preserves the structure of the canonical Poisson 
brackets 

{N,U N } = l = {u,U v }, (50) 
{X",p„} (51) 

as expected. 

In terms of the new phase space variables, the first- and second-class constraints 



(36)-(39) become 
NUm = 0, 



and 



A' 



/r 



Pt + /3 — cosh n„ + (p r — 2arU u ) sinh IL, + 2ar cosh n„ 



S 2 



0. 



N 



p t + (3 — ) sinh n,, + (p r — 2arU u ) coshll^ 



(52) 
(53) 

(54) 
(55) 



respectively. Note that in the second-class constraint (55) we split the momentum 



conjugated to r according to relations (19) and (49). As customary, the second-class 



constraints (54) and (55) may be taken as algebraic identities after implementing the 



Dirac bracket [13 [20] ■ Furthermore, these second-class identities will become auspicious 
at the quantum level since they enclose important operator identities. 

The constraint T\ is simply associated to the gauge transformations Nd^ — Hn9u n 



which only acts on the A^IlAr-plane. As for the constraint T% in equation (53), we can 



further transform it by expressing the hyperbolic functions in terms of the phase space 
variables as follows, 



cosh n. 



sinh n,. 



Pt + P- 



r 2( 2a7 )i/2 

(p r - 2arU. v ) 
r 2 {2ai) x l 2 



(56) 
(57) 



where we have substituted the second gauge condition (44) into the second-class 

/VV n f , ^ 



constraint (|55|). Thus, T i is transformed into 

•^"2 o /n \ i /o (P r 



! (2a 7 ) 1 /2 



2arU u y - [p t + 



2ar p t + 



0. 



(5f 
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and we have arrived to an expression quadratic in the physical momenta for the canonical 
Hamiltonian H Q , which is identified with the constraint Ti when the second-class 
constraint S\ is considered. 

In order to remark the relevance of the new canonical variables, we observe that the 



internal energy Q given by expression (17) can be written as Q = 2orcosh U u + (3q 2 /r 



which at first-order approximation matches the potential energy in [1] in the so-called 
low velocity regime. However, it is easy to show that an additional kinetic term 
emerges from this expression by considering up to second-order in the expansion 



2 M (r) 



of coshllj,, where /i (r) = (4ar) -1 . Note that this mass-like term is a function 
of r. In addition, if we consider the first-class constraint Ti y we can obtain the 
expression Q = ±a/ (p r — 2arU u ) 2 + m 2 (r) + V re \(r), where the mass-like term is given 

by m 2 (r] 



2^,2 



a r 



and the relativistic potential is given by V re i(r) := /3q /r + ar. 



7. Quantum approach 

In this section we study the canonical quantization of our system. To this end, we 
emphasize the totally dissimilar nature which first- and second-class constraints play in 
the quantum theory. We start with the conventional way by promoting the classical 
constraints into operators. However, by implementation of the Dirac bracket in the 
classical counterpart we are able to eliminate second-class constraints off the theory by 
converting them into strong identities. At the quantum level this is mirrored by defining 
the quantum commutator of two quantum operators, A and B, as 

[A, B] :=i{XB}*, 

where {•, •}* stands for the Dirac bracket. Thus, with this prescription the operators 
corresponding to second-class constraints are also enforced as operator identities [IH [20] . 
For the system in question, this yields the quantum operator expressions 

S 1 :=u = 0, (59) 

7q 2 \ 



N 



p t + (3— J sinh Ii u + (p r — 2arIL y )cosh n,, 







(60) 



which, in particular, tell us the character of the quantum operators z>, coshll^ and 
sinh H",, (see equations (56) and (57)). Also, we will represent the radial operator as 
p r := —i(d/dr)r since then the operator p 2 . = —{d 2 /dr 2 + (2/r)d/dr) will be Hermitian 
in the inner product of states in a conventional Hilbert space, namely an L 2 -space. For 
the rest of the variables, we choose to work on the "position" representation, where 
we consider the position operators by multiplication and their associated momenta 
operators by —i times the corresponding derivative operator. 

Next, we will adopt as our quantum first-class constraints the operators 

p l ■= - iN^r , (61) 



T, 



N 



dN 



p r - 2arn i/ ) - [pt + 



/V 



— 2ar ( p t + 



/V 



(62) 
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Note that the N factor in Ti is necessary in order to maintain at the quantum level the 
classical algebraic structure between the two first-class constraints. Also note, that for 
simplicity, we have chosen a trivial factor ordering which allowed us to get rid of the 
denominator in equation (p^p] First, we explore the quantum equations emerging by 
considering the physical states \1/ of the theory as those defined by naive Dirac conditions 



T 2^ 



0, 
0. 



(63) 
(64) 



Equation (63) simply tells us that our physical states \I/ are not explicitly depending on 



the phase space variable N. As classically, equation (64) is the most interesting for us 



since it is related to the Hamiltonian operator Hq, hence resulting in a Schrodinger-like 
equation. In order to find solutions to this equation, we notice first that the second- 
class constraint S\ tell us that the variable v is fbced to zero, thus getting rid of any 
possible dependence on this variable in which in turn leads to the conclusion that the 
action of the operator il„ on the \1/ states vanishes automatically. Further, we see that 
the t-dependence can be solved by assuming in aereement with 

the classical definition for Q. Finally, we notice that ip(r) must satisfy the differential 
equation 



d 2 2 d 
dr 2 r dr 



V{r) 



ij){r) = 



where the quantum potential V (r) is given by 
V (r) = 



P 2 q 4 2PQq 2 
— 7T- H h 2a\lr 



n 2 - 2a/3q 2 



(65) 



(66) 



The shape of this potential is depicted in FigureJT} The real zeroes of V (r) are located at 

if the condition 



n > ' 2 



H. 1 
4o 1 



1 



jp-j and r 3 :=£^l + ^l-^ 

Q 2 > 8af3q 2 holds. We will work under this condition. It follows from the preceding 
expressions for the zeroes of the potential that if Q increases, n and r 2 come to zero 
and, in consequence, the local maximum decreases. In opposition, as Q increases the 
zero located at r 3 grows. The potential tends to infinity as r grows whereas it becomes 

+ Indeed, if we consider a different factor ordering for the quantum counterparts of first-class 



constraints, (52 1 and (58), we have as a result that the wave function is a homogeneous function of 



degree minus one half in the iV-variable, and also, the potential depending on the r-variable (described 



below in (66 1 ) includes an extra inverse square term. Important to mention is the fact that the common 



denominator in ( 58 ) is cancelled out in the low velocity regime of [3] , and hence the factor ordering they 



have chosen results simpler in that case. From a more extensive point of view, the operator ordering 
problem can be tackled by following the arguments developed in quantum cosmology as done in |21| 
and [22 , where the operator ordering is related to the avoidance of singularities and to the Hermiticity 
of the Hamiltonian operator. In fact, among the most popular choices (not unique) for a suitable factor 
ordering in quantum cosmology, is to consider expressions with derivatives in the form of a Laplace- 
Beltrami operator. At a fundamental level in our quantum description there is not obvious election for 
a factor ordering. Our choice for the operator p r results Hermitian and in consequence the implicit 



operator ordering ambiguities are fbced in the expression (62) (for example, compare with Eq. (A.l) 
in |21j with p = 2 in flat minisuperspace). 
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Figure 1. Effcctive quantum potential of a charged rigid bubble for large values 
of g and when g approaches to zero. For tiris graph we considered the energy value 
e = 0.1144817. 



singular at minus infinity as r is approaching to zero. Note that this effective potential 
has a linearly raising dependence of the bubble radius which is different in comparison 
with the original Dirac model which possesses a quartic dependence in the bubble radius 
added with the first, the second and the fourth term of the potential (66). The first 
term in the potential (66) corresponds to a centrifugal-like potential whereas the second 
term is related to the usual Coulomb interaction and, the latter term may be conceived 
as a shift in the bubble energy levels by the quantity 2a(3q 2 . 

The issue of computing explicit solutions of the Eq. (65) results very complicated, 
instead we bring into a play a numerical technique based in a Fortran code in order 
to compute the energy spectrum. For this reason, it is convenient to choose the 
dimensionless quantities, g := r/r and e = r Q (see for example [23]) where we have 



considered the results of the Subsection 4^ Now, by a suitable alternative function 
0(r) := ip(r)/r, the radial equation (65) becomes 

d 2 



d g 2 



where 



V{g) 



+ V(q))<Kq) = 0. 



X 2 e\ e\ 



4g< 



+ h 



e 2 + 



4 



(67) 



(68) 



By employing the fourth-order Runge-Kutta method altogether with the necessary 
boundary conditions 0(0) = and (j)'(0) = A, the latter being an arbitrary value which 
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Figure 2. The first normalized eigenfunctions of a charged rigid membrane. 



is not relevant for the final result, we are able to obtain the solutions of the Eq. (67) 



Now, by using the shooting method, both the eigenfunctions and eigenvalues of the 



Schrodinger-like equation (67) can be computed. Some of the excited states eigenvalues 



of the charged bubble are shown in Table [T] Similarly, in Figure [2] we plot the first 
normalized eigenfunctions <f) n (g). At this point, we are able to compare our results with 
those developed by Onder and Tucker in j I]. We found that for n = 1, Qi = 15.68 m e 
and for n = 133, ^133 = 206.768 m e (the muon mass), unlike the values obtained by OT 
where, for instance, they found for the ground state n = 1, Oi ~ 13 m e and for n = 190, 
f2i9o — 200 m e . This fact is a merely consequence of the several aspects of the quantum 
approach developed by OT where a semiclassical quantization, in a specific gauge, was 
adopted. Hence, the old idea pursued by Dirac regarding that the first excited state 
could be considered as a muon is not realized for this rigid bubble model. 



8. Concluding remarks 

In this paper we reviewed the classical and quantum aspects of a rigid charged 
conducting membrane. We considered our model as a dynamical bubble in the presence 
of an electromagnetic field with the non-electromagnetic forces described by a term 
depending linearly on the extrinsic curvature of the worldvolume. As discussed before, 
though such correction term involves second-order derivatives, the ensuing equations of 
motion are of second-order in the field variables, as expected. We would like to remark 
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Table 1. Eigenvalues of the Schrodinger-like equation (67 1 and the excited energies 
of the charged bubble 



n 




Q n /m e 


1 


0.1144817 


15.6839929 


2 


0.1738982 


23.8240616 


3 


0.2177701 


29.8345146 


4 


0.2541984 


34.8251903 


5 


0.2860303 


39.1861634 


G 


0.3146605 


43.1085008 


7 


0.3408957 


46.7027122 


8 


0.3652517 


50.0394897 


9 


0.3880928 


53.1686423 


10 


0.4096428 


56.1210636 



that a generic minimal coupling term between the charge and the field has been 
considered in this work [11], and it is equivalent to the Dirac original proposal that 
introduces a coupling by means of a boundary condition which is consistent only in a 
particular gauge p]. 

Our analysis took into consideration the routinely neglected boundary term which 
involves second-order derivatives of the configuration variables. Thus, we evoked the 
Ostrogradski-Hamiltonian formalism in order to investigate the constrained structure 
of the model. Indeed, a canonical transformation allowed us to acquire control over 
our set of constraints.. As expected, we are ushered to a Hamiltonian constraint which 
do not depend on the momenta associated to the higher-order variables, hence the 
structures coming from the boundary term do not play a relevant part on the dynamics 
of the system, but nonetheless they play the important role of being a bridge to obtain 
quadratic expressions in the physical momenta, useful for ge to the quantum 

theory. 

At the quantum level, we considered the canonical formalism of Dirac in order 
to overcome certain issues found in previous attempts by several other authors. In 
particular, we obtained a Schrodinger-like equation where the explicit effective quantum 
potential is identified, and the global behavior of the potential is discussed. As a 
byproduct, we have made a reappraisal of the excited states eigenvalues of this type of 
charged membrane which demonstrates the significant value of prediction of the scheme 
followed here. 

More comments are in order. The Dirac model presents signals of instability [25] . 
Perhaps, this shortcoming could be overcome if we take into account a linear extrinsic 
curvature correction term in the model. Other proposals to eliminate the stability 
problem of the original Dirac model have been approached by considering a spinning 
bubble in the way studied in [26] or by the introduction of worldvolume fields pondering 
the spin of the electron [27]. Certainly, if a more complete model which includes the 
tension of the membrane is considered, we will not necessarily bring an improvement on 



Quantum charged rigid membrane 



16 



the dynamical behavior since it will engender some singularities or discontinuities in the 
potential which can not be ruled out from the outset. This caveat must be scrutinized 
cautiously in the dynamical analysis and provides a good motivation for subsequent 
studies. 

Finally, the model described here, results interesting in its own right since it has 
all the hallmarks of cosmological brane models. Among the brane cosmology concerns, 
we can mainly point two of them: their gravity and their dynamics. In this context, 
making use of the scheme here developed, it will be interesting to study the cosmology 
associated to the Dvali-Gabadadze-Porrati model, now taking into account a term 
depending linearly on the extrinsic curvature with the purpose to remove the pathologies 
and instabilities of such model [28]. An enormous advantage will be that the resultant 
effective model will retain a linear dependence on the acceleration of the field variables 
which open the possibility to apply the same quantum strategy followed here. 
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